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Abstract 

The line of phase equilibrium between the hot and cold phases of a plasma is derived from 
the intensity of thermal radiation as a function of the plasma density and temperature. The last 
function is obtained with the help of the condition for emission relating the wavelength of radiation 
and the plasma density. 
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It was shown recently [1] that a transport barrier in toroidal plasma, i.e. the region of 
reduced energy transport across magnetic surfaces [2-4] , is the boundary region between the 
hot and cold phases of the plasma. The transport properties of the hot phase essentially 
differ from those of the cold phase due to the intense interaction of the plasma with thermal 
radiation above the inversion temperature [5] . The boundary region between the hot and cold 
phases corresponds to the minimum values of specific resistivity and transport coefficients 
of the plasma. Other specific properties of a transport barrier such as low magnetic shear, 
the generation of radial electric fields, and the poloidal rotation are the consequences of the 
minimum resistivity and transport coefficients of a plasma in this region. 

The specific resistivity of the cold phase in the vicinity of the phase transition to the hot 
phase is very small, so the phase transition is nearly superconductive. The additional current 
flows in the cold phase nearby the boundary with the hot phase and does not heat the hot 
phase, so the additional neutral-beam or radio-frequency heating is required to obtain higher 
temperatures in the hot phase [3,4]. 

The critical value T e c of the electron temperature, corresponding to the transition from 
the cold to the hot phase, for relatively low plasma densities achievable in tokamaks is 
independent of the density and has an order of magnitude of the inversion temperature 
T = 2keV [5]. Studies of sawtooth oscillations by von Goeler, Stodiek, and Southoff [6] 
suggest that the transition temperature is approximately T 6iC = 0.8keV. In this paper, we 
expand the line of phase equilibrium between the hot and cold phases of a plasma in the 
region of higher densities where the dependence of the transition temperature on the plasma 
density is essential. 

Earlier, the energy distribution of atoms in the field of thermal blackbody radiation was 
obtained [7] in the form 

jV/jV = a a u 2 / (2ttc 2 ) (exp (hw/kT) - 1) , (1) 

where Nq is the population of the ground state Eq, N is the population of the energy level 
E = E + ftu), a a is the absorption cross-section, h is the Planck constant, and T is the 
radiation temperature. 

This distribution is valid in the range hw/kT ^ 1, since in the limit hw/kT — ► the line 
width is going to infinity, that indicates the violation of the one-particle approximation used 
in [7]. 
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The function has a maximum at hu m = l.GkT. When the temperature exceeds the 
critical value of To = 2 x 10 7 K (the inversion temperature), the population of the energy 
level E exceeds the population of the ground state E . Since the function (JTJ) is increasing 
in the range u < u m , the inversion of the energy level population is produced also in 
some vicinity of uo m (below uj m ). This suggests the maser amplification of thermal radio 
emission in continuum by a hot plasma with the temperature exceeding the critical value 
To. Maser amplification in compact radio sources was assumed earlier in [8] based on the 
high brightness temperatures of active galactic nuclei. Since a hot plasma in an accretion 
disk is concentrated nearby the central energy source, maser amplification is characteristic 
for compact radio sources. 

It is clear that, when the temperature of a plasma is below T , the radio flux is very small, 
and when the temperature exceeds To, radio emission is on. This an on-off cycle is detected 
in the radio pulsar PSR B1259-63 [9]. Similar is an on-off cycle in X-ray pulsars, e.g., the 
35-day cycle in Her X-l. It implies that X-ray emission from X-ray pulsars is produced by 
the laser amplification in continuum which is quite analogues to maser amplification at radio 
wavelengths. 

It was shown also [7] that thermal emission has a stimulated character. According to this 
conception thermal emission from non-uniform gas is produced by an ensemble of individual 
emitters. Each of these emitters is an elementary resonator the size of which has an order 
of magnitude of mean free path I of photons 

l = — (2) 
no- 
where n is the number density of particles and a is the absorption cross-section. 
The emission of each elementary resonator is coherent, with the wavelength 



A = I, (3) 

and thermal emission of gaseous layer is incoherent sum of radiation produced by individual 
emitters. 

An elementary resonator emits in the direction opposite to the direction of the density 
gradient. The wall of the resonator corresponding to the lower density is half-transparent 
due to the decrease of absorption with the decreasing gas density. 
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The condition (J3J) implies that the radiation with the wavelength A is produced by the 
gaseous layer with the definite number density of particles n . 

The condition (JHJ) is consistent with the experimental results by Looney and Brown on 
the excitation of plasma waves by electron beam [10]. 

The intensity of thermal blackbody radiation is given by the Planck formula 

B v (T) = 2hu 3 /c 2 (exp (his/T) - l)" 1 , (4) 

where v = c/X is the frequency of radiation, c is the speed of light, h is the Planck constant, 
T is the temperature, the Boltzmann constant being assumed to be included in the definition 
of the temperature. Using the condition for emission (j3J), we obtain the intensity of thermal 
blackbody radiation as a function of the density, n, and the temperature, T, in the form 

B (n, T) = (2T 3 //i 2 c 2 ) {an/Tf {exp {an/T) - , (5) 

where a = ^hca a , anda a is the absorption cross-section. 

Since the function / (x) = x 3 / {e x — 1) has a maximum at x m ~ 3, the value of the 
density, n c , corresponding to the maximum value of the intensity of thermal radiation, is 
given by the formula 

n c « 3T/a = 12T/ {hca a ) . (6) 

A hot plasma with the temperature T e ^ T = 2keV is intensely interacting with the field 
of thermal radiation. At temperatures T ^ To the stimulated radiation processes dominate 
this interaction [5]. Thermal radiation induces radiative transitions in the system of electron 
and ion which corresponds to the transition of electron from the free to the bounded state. 

Thus, in a hot plasma interacting with thermal radiation, the bounded states of electrons 
and ions restore, leading to the change of collision properties of a hot plasma. In this case, 
the electron-ion collision cross-section has an order of magnitude of the atomic cross-section, 
a = l(T 15 cm 2 . 

In the Wien region n > n c (T) the intensity of thermal radiation exponentially decreases 
with the increase of the density, so the interaction between a plasma and thermal radiation 
is weak, and hence this region corresponds to the cold phase of a plasma. The hot phase of a 
plasma is realized in the Rayleigh- Jeans region n < n c (T) (assuming that the temperature 
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exceeds the inversion temperature), where the intensity of thermal radiation comparatively 
weakly depends on the plasma density. 

As a result, the line of phase equilibrium between the hot and cold phases of a plasma is 
described by the following approximate expressions: 

T c = T C)0 , n < n Cj0 , (7) 

n c ^12T/(hca a ),n>n cfi , (8) 

where n Cj0 = 0.8 x 10 17 cm~ 3 for T Cj0 = 0.8keV, according to the equation (6). 

Consider now the phase equilibrium between the hot and cold phases of a plasma. The 
pressure in the cold phase is given by the formula 

P c = 2mT, (9) 

where n« = n e are the ion and electron densities respectively. The pressure in the hot phase 
is 

P h = n T + 2(n h -n )T, (10) 

where n is the density of neutral atoms, and is the total density of ions and neutral 
atoms. 

If no ~ rih, then from the condition of equilibrium, P c = Ph, we find that ~ n c , i.e. the 
density of the hot phase is approximately twice of the density of the cold phase. It means 
that there is the discontinuity of density at the boundary between the hot and cold phases, 
analogous to the first order phase transitions. 

There is also the discontinuity of resistivity somewhat analogous to the superconduction 
phase transitions. The specific resistivity of the hot phase, rjh, is about three orders of 
magnitude larger than that of the cold phase, i] c . As much as large is the discontinuity of 
thermal conductivity across magnetic surfaces, so the cold phase is a good thermal protector 
for the hot phase. 

The discontinuity of diffusion coefficients with respect to the diffusion across magnetic 
surfaces leads to the generation of radial electric fields in the region of a transport barrier 
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due to the differences in diffusion coefficients for ions and electrons . The drift of electrons 
and ions in these electric fields produces the poloidal rotation of the plasma. 

To summerise, the critical temperature of the phase transition between the cold and 
hot phases of a plasma is independent of the plasma density for low densities and increases 
proportionally to the plasma density for higher densities exceeding some critical value. There 
are discontinuities of density, resistivity, and transport coefficients at the point of phase 
transition. The high-temperature phase transition in a plasma is experimentally observed 
transport barrier in toroidal plasma. 
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